Ferroelectricity in spiral magnets 
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It was recently observed that materials showing most striking multiferroic phenomena are frus- 
trated spin-density-wave magnets. We present a simple phenomenological theory, which describes 
the orientation of the induced electric polarization for various incommensurate magnetic states, its 
dependence on temperature and magnetic field, and anomalies of dielectric susceptibility at magnetic 
transitions. We show that electric polarization can be induced at domain walls and that magnetic 
vortices carry electric charge. 
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The recent revival of interest in materials showing 
strong interplay between magnetism and ferroelectricity 
led to the discovery of a new class of systems, which can 
be called ferroelectric magnets. While in 'old' multifer- 
roics, such as BiFe03, spontaneous electric polarization 
appears at a much higher temperature than magnetism, 
in RMn0 3 (R = Tb, Dy, Gd) 0,110, Q, RMn 2 5 (R 
= Tb, Ho, Dy)H@] and Ni 3 V 2 8 0] ferroelectricity is 
only observed in magnetically ordered states. In mag- 
netic fields these materials show reversals and sudden 
flops of electric polarization vector an d an excep- 

tionally strong enhancement of dielectric constant (the 
giant magnetocapacitance effect) 3]. These remarkable 
multiferroic phenomena, essential for control of dielectric 
properties by magnetism, follow from the fact that that 
ferroelectricity in these materials is induced by magnetic 
ordering. 

An intriguing feature of the new class of multiferroics is 
an intimate link between the spontaneously induced po- 
larization and magnetic frustration. In RMn03 a large 
distortion of the cubic perovskite lattice gives rise to com- 
peting ferro- and antifcrromagnetic interactions between 
the Mn spins [l| , while in the Kagome staircase material 
M3V2O8 frustration originates from the lattice geometry 
0- The competing exchange interactions stabilize spin- 
density-wave (SDW) states with a periodically varying 
magnetization 

The relation between the ferroelectricity and incom- 
mensurate magnetism is widely used as an empirical 
guiding principle in the search for new multiferroic ma- 
terials. In this paper we discuss a phenomenological de- 
scription of ferroelectric magnets, based on general sym- 
metry arguments. We formulate a simple continuum 
model, which clarifies the relation between the induced 
electric polarization and magnetic structure, describes 
anomalies of dielectric constant at magnetic transitions 
and qualitatively explains complex magnetic field behav- 
iors found in these materials. We extend this analysis 



to domain walls and vortices. Our phenomenological ap- 
proach is complimentary to the recent discussions of mi- 
croscopic mechanisms of ferroelectricity in magnets 0,0- 
Induced polarization: Incommensurate SDW states are 
largely insensitive to details of crystal structure and can 
be described by a continuum field theory of the Ginzburg- 
Landau type. The form of the coupling of electric polar- 
ization P to magnetization M can be found using general 
symmetry arguments. The invariance upon the time re- 
versal, t — > —t, which transforms P — » P and M — > — M, 
requires the lowest-order coupling to be quadratic in 
M. The symmetry with respect to the spatial inversion, 
x — » x, upon which P — * P and M — > M, is respected 
when the coupling of a uniform polarization to an inho- 
mogeneous magnetization is linear in P and contains one 
gradient of M. Omitting vector indices the coupling term 
in thermodynamic potential can be written in the form 



4> em (P,M) oc PMdM. 



(1) 



The terms linear in gradient (Lifshitz invariants) are al- 
lowed in systems with broken inversion symmetry, such as 
noncentrosymmetric crystals, where they can give rise to 
periodic spatial modulations of magnetization [loj . Such 
an incommensurate SDW state is observed in the ferro- 
electric BiFeOs, where the inversion symmetry is spon- 
taneously broken by electric polarization. The helix with 
the long period 620A in BiFeOs results from the coupling 
between the electric polarization P and the G-type AFM 
order, which has the form Eq.Q with the uniform mag- 
netization M replaced by the slowly varying Neel vector 
L The small value of the wave vector is a conse- 

quence of the relativistic nature of the coupling 0, 0] . 

This reasoning can be turned around to explain electric 
polarization in frustrated magnets, in which an incom- 
mensurate magnetic ordering results from competing ex- 
change interactions and the SDW wave vector is, in gen- 
eral, not small. When the SDW order of a proper kind 
sets in, the coupling Eq.QJ induces a uniform electric 
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polarization that breaks the inversion symmetry. The 
weakness of the coupling translates in this case to rela- 
tively low values of the induced polarization. 

This mechanism does not require a special kind of crys- 
tal lattice. In the simplest case of cubic symmetry the 
coupling term has the form 

$ em (P, M) = 7 P • [M (V • M) - (M • V) M + . . .] . (2) 

The omitted terms can be written as the the total deriva- 
tive, V/(M), and do not contribute to the uniform po- 
larization. Assuming that in absence of magnetism the 
system shows no instability towards ferroelectricity, we 
only keep the quadratic term in the 'electric part' of the 
thermodynamic potential, 3> e (P) = where \ e is the 
dielectric susceptibility in absence of magnetism. The 
variation of $ e + <i> em with respect to P then gives 

P = IXe [(M ■ V) M - M (V • M)] . (3) 
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FIG. 1: The sinusoidal SDW [panel (a)] does not induce a uni- 
form electric polarization, while for a helicoidal SDW [panel 
(b)] the electric polarization vector P lies in the spin rotation 
plane (with the normal vector es) and is transverse to the 
wave vector Q of the helix. 

Consider now a SDW state with the wave vector Q, 

M = A/iei cos (j) + M2G2 sin <f> + M363, (4) 

where <f> = Qx and the unit vectors e^, i = 1,2,3 form an 
orthogonal basis. If only Mi or M 2 is nonzero, this state 
is the sinusoidal wave (see Fig. [TJi) , while for M\ , M 2 ^ 
it is an (elliptical) helix with the spin rotation axis e3. 
If also M3 ^ 0, then the helix is conical. Using Eq.©, 
we find that the average polarization is transverse both 
to the spin rotation axis e3 and the wave vector Q of the 
helix (see Fig. QJd), and is independent of M3: 

P = I j d 3 xP = -fXeMxMi [e 3 x Q] . (5) 

For the sinusoidal SDW the induced polarization is 0: 
such an ordering does not break the inversion symmetry 



at the sites where the magnetization reaches maximum or 
minimum (e.g., the points on the dashed line in Fig.^i) 
and, therefore, it cannot induce an average electric po- 
larization. Equation (J5J also holds for the orthorombic 
crystal symmetry, provided that e3 and Q are parallel to 
crystal axes. 

This explains why the transition at T$ — 41K to the si- 
nusoidal SDW state in TbMnO^[2lH^| does not give rise 
to ferroelectricity. The polarization is only induced below 
the so-called lock-in transition at Th = 28K, when the 
sinusoidal SDW is replaced by the helix with the propa- 
gation vector Q parallel to the b axis and the Mn spins 
rotating in the be plane (e 3 || a) [Hi}, so that according 
to Eq.JSJ the polarization is induced along the c axis in 
agreement with experiment. Similarly, the polarization 
is absent in the high-temperature incommensurate phase 
of Ni 3 V 2 8 , which is the sinusoidal SDW state 7J. The 
helix with Q || a and the spin rotation axis e3 || c, which 
appears in the low-temperature incommensurate phase, 
induces polarization along the b axis. 

We note that the spiral ordering can be considered as 
a particular case of a magnetic ordering with two non- 
collinear SDWs with equal wave vectors, 

M = Miei cos (Q-x + fa) + M 2 e 2 cos (Q • x + fa) , 

such as the one recently found in RMn 2 05 [3E3- Al- 
though a single sinusoidal SDW does not induce polar- 
ization, the interference between the two SDWs gives 

P = jXeMxMi sin (fa - fa) [Q x [e x x e 2 ]] . (6) 

Magnetic textures: Electric polarization can also be in- 
duced in nonfrustrated magnets near magnetic defects, 
e.g., domain walls or in inhomogeneous ground states 
stabilized by magnetostatic interactions, e.g., vortices in 
nanodiscs 16]. In many cases the integrated quantities, 
such as the total polarization per unit area of a domain 
wall, only depend on topology of spin textures. For walls 
with collinear spins and the ones with spins rotating 
around the axis normal to the wall, the total polarization 
is 0, while the wall, in which spins rotate around an axis 
parallel to the wall, M = M [cos faxi) ei + sin</>(a;i) e2], 
has the total polarization that depends on the total ro- 
tation angle 

/+00 
dxiP = 2 7Xe M 2 e 2 [fa+oo) ~ fa-oo)\ . (7) 
-00 

Using Eq.© one can also show that the vortex 

M = M [cos(n</> + (/>o)ei + sin(n</> + (/>o)e 2 ] , 

where (f> — arctan ^ and fa is an arbitrary phase, has 
quantized electric charge located at the vortex core: q n = 
nqi, where qi — AirjXeM 2 and n is the winding number 
of the vortex. An applied electric field will move magnetic 
vortices and anti-vortices in opposite directions. 
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Sinusoidal-helicoidal transition: We now turn to the 
phase diagram of ferroelectric SDW magnets. Using the 
values of the induced polarization (10 2 — 10 3 ^C m~ 2 ) 
and magnetic transition temperatures (5 — 40K) for these 
materials, we find that the energy gain related to the 
induced polarization is small compared to the magnetic 
energy gain. Therefore, the temperature and magnetic 
field dependence of the polarization merely reflects the 
changes in magnetic ordering, which can be described 
using the Ginzburg-Landau thermodynamic potential 



-M 



-M 



t—x,y,z 



dx 2 



Q ) M. 

(8) 

In what follows we assume that a x < a y < a z (easy 
axis along the x direction) and first neglect higher-order 
anisotropies. The last term in Eq.JHJ) favors a periodic 
SDW ordering with the wave vector Q along the x axis. 

A down-shift of the ferroelectric transition with respect 
to the magnetic one, found in all magnetic ferroelectrics, 
is a consequence of magnetic anisotropy. While for an 
isotropic system the ground state is a helix with a con- 
stant M, an anisotropic system first undergoes a transi- 
tion to the sinusoidal SDW state with M along the easy 
axis, M = M x x cos Q x, at temperature T$ : a x {T$) = 0. 
As temperature is lowered and the amplitude of the order 
parameter grows, the system undergoes a second transi- 
tion at some Tjj < T$ to the elliptical helix state, M = 
M x x cosQx + M y y smQx, provided that the anisotropy 
parameter A = not too large. When the two 

transitions occur at close temperatures, the higher har- 
monics in the SDW state are small and the helix appears 
at a y = Or/3. For a x (T) — a(T — T$) we then obtain 

T H = T S -^. (9) 

The average electric polarization only appears in the 
helicoidal state and for spins rotating in the xy plane and 
Q || x it is parallel to the y axis: 



Py = otlXeQV (Th - T) (T S + A/(2a) - T). (10) 

Note that since P y cx M x M y [see Eq.ljSJ)], it has the 
square root anomaly at the ferroelectric transition, even 
though it is not a primary order parameter. Further- 
more, as in proper ferroelectrics, the dielectric constant 
e yy diverges at Th and obeys 'the 1/2-law' [r^ : 



T-Th ■ 



-yy 



for T > T H , 



r , for T<T H , 



(11) 



2(T H -T) : 

where A = 6A(7y e Q) 2 /(a&). Although P y vanishes in 
the sinusoidal SDW state, the magnetic contribution to 
e yy is nonzero up to T = Ts' 



(2 1X eQ? (T s - T) 
b (T - T H ) ' 



for T H <T< T s . (12) 
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FIG. 2: Magnetic field behavior of electric polarization for 
the model Eq.®. In zero field spins rotate in the xy plane 
and P || y. Magnetic field in the x direction suppresses the 
polarization, while H y orients P in the z direction. 



Behavior in magnetic field: The salient feature of fer- 
roelectric magnets, important for technological applica- 
tions, is the strong sensitivity of their dielectric properties 
to magnetic field, which can suppress electric polarization 
or change its direction 0, Q . We first discuss polariza- 
tion flops in the model Eq.©. In weak fields spins rotate 
in the easy xy plane, so that the spin rotation axis e3 of 
the helix is parallel to the 'hard' z axis. For Q || x elec- 
tric polarization is oriented along the y axis. In strong 
magnetic fields spins form a conical helix with e3 || H, 
e.g., H x will force the spins to rotate in the zy plane (see 
Fig. 130- Such a spin flop will suppress electric polar- 
ization, since for the zy-helix e3 || Q and according to 
Eq.© P = 0. On the other hand, magnetic field in the 
y direction favors the rotation of spins in the xz plane, 
in which case P || z (see Fig.|2p). 

The magnetic field behavior observed in orthorombic 
manganites is somewhat more involved. If we identify 
the x, y, and z axes used in this paper with, respectively, 
the 6, c, and a axes of the Pbnm crystal structure of 
TbMnOa , then magnetic field applied in the x and z 
directions changes the direction of the electric polariza- 
tion of TbMn03 from y to z. According to Eq.JSJl, this 
corresponds to the change of the rotation plane from xy 
to xz. It is the flop shown in Fig. |2p, but induced by 
magnetic fields with 'wrong' orientations. 

This unusual behavior is most likely related to the flops 
of the strongly anisotropic rare earth spins, coupled to 
Mn spins [3.ll3| . It can be described phenomenologically 
by adding the higher-order anisotropies to Eq.©, e.g., 



A<r V (M) = l, lll{ \[ ir iM y ) 2 + b' xy (M x ) 2 (^ 



+ b yz {M y f{M z f 



(13) 



For positive coefficients the first two terms suppress the 
rotation in the xy plane, when magnetic field is applied 
in the x or y direction, while the last term suppresses 
the rotation in the yz plane. This gives rise to phase 
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diagrams, shown in Fig. which are similar to the ones 
found for TbMnOa jj| . The nonlinear terms also result in 
a magnetic field dependence of the helicoidal transition 
temperature, which together with the divergency of di- 
electric constant at Tjj makes e strongly field-dependent 
(the giant magnetocapacitance effect |3|). 
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of P(H X ), using the Ginzburg-Landau expansion similar 
to Eq.©, is shown in Fig. 0] 




FIG. 4: Predicted magnetic field dependence of the electric 
polarization \P V \ = |P Z | of CaFeOa for H || x. 

In conclusion, we used simple symmetry arguments to 
explain ferroelectric properties and thermodynamics of 
spiral magnets. Taking into account the complexity of 
exchange interactions and spin orders in these materials, 
this phenomenological approach works surprisingly well. 
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port by the DFG (Merkator fellowship) and MSC P ' US pro- 
gram is gratefully acknowledged. 
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FIG. 3: Typical phase diagrams of the model Eq.JSJ with 
nonlinear terms Eg. 11311 included for H || x, y,z [respectively, 
panels (a),(b), and (c)]. Here X-SDW (green) denotes the 
sinusoidal SDW state with spins along the x axis, while XY- 
SDW (pink) and XZ-SDW (blue) denote the spiral states with 
spins rotating, respectively, in the xy and xz planes. 

Importantly, magnetic fields required to induce spin 
flops are of the order of magnetic anisotropies that can be 
relatively small for transition metal ions with filled and 
half-filled t 2g shells, e.g., Cu 2+ , Ni 2+ , Fe 3+ , and Mn 3+ . 
Therefore, electric polarization can be flopped by modest 
magnetic fields even in spiral magnets with high ordering 
temperatures, which may be interesting for applications. 
Consider, e.g., a gedanken experiment on the helimag- 
net CaFe0 3 with T N = 115K In zero field both 

the wave vector and the spin rotation axis of the helix 
are parallel to the body diagonal, Q,e3 || [1,1,1], and 
no electric polarization is expected. Magnetic field of 
the order of magnetic anisotropies can flop the orienta- 
tion 63. On the other hand, to change the wave vector 
Q would require fields of the order of the antiferromag- 
netic superexchange between neighboring iron spins |l9| , 
which are much stronger. Thus magnetic field along one 
of the crystal axes, e.g. H = H x [1,0,0], will induce 
P = P[0, —1, 1]. The result of the numerical calculation 
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